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Abstract 

The interplay of spectroscopy, scaling laws and critical indices is 
studied in strongly coupled quenched QED. Interpreted as a model of 
technicolor having strong interactions at short distances, we predict 
the techni-nieson mass spectrum in a simplified model of a dynami- 
cally generated top quark mass Mf. Our results support the strict 
inequality that the techni-sigma mass Ma- is less than twice the dy- 
namical quark mass Mf , and confirm that the techni-pion is a Nambu- 
Goldstone boson. The level ordering = M^r < M^ < 2Mf < Mp < 
Mai is found. An equation of state, and scaling laws are derived for 
the techni-meson masses by exploiting correlation length scaling. The 
resulting universality relations are confirmed by simulations on 16^, 
32 X 16^ and 32^ lattices. The anomalous dimension t] is measured to 
be approximatively 0.50 in good agreement with past lattice simula- 
tions and hyperscaling relations, as well as with the analytic solution 
of the quenched, planar gauged Nambu-Jona Lasinio model solved by 
continuum Schwinger-Dyson equation techniques. 



1 Introduction 

Quenched, strongly coupled QED is an interesting, relevant model field 
theory |l| , . Since the photon sector of the theory is free, the model is 
amenable to analytic solution and a wealth of interesting dynamics has been 
found in its light fermion sector. Chiral symmetry breaking, a dynamical 
realization of the Goldstone mechanism, a rich renormalization group struc- 
ture, large anomalous dimensions, super-critical field phenomena, etc. have 
been discussed and illustrated in this context I]!], 0, |§[, |0], 0, 0- 

Many of the continuum predictions have been confirmed by lattice simu- 
lations which have stressed the theory's chiral Equation of State (EOS) [0]. 
The central purpose of the present paper is to study the theory's mass spec- 
trum, and relate the spectrum to other fundamental properties of the theory, 
such as its critical indices (anomalous dimensions). Since the model is a 
simplified version of Technicolor dynamics f^ and has been advocated as the 
dynamical basis of the mass of the top quark, our results should be of interest 
to high energy model builders looking beyond the Standard Model toward the 
physics underlying the Higgs mechanism |^. Using lattice techniques bor- 
rowed from lattice QCD spectrum calculations we shall calculate the masses 
of the techni-fermion(/), pion (vr), rho (p), sigma (a) and al (al) [0. In 
the strong coupling phase where chiral symmetry is broken dynamically, we 
confirm that the pion is massless and satisfies the constraints of PCAC The 
other states are massive with the ordering Mf < Ma- < Mp < Mai. In fact. 
Ma- < 2Mf as expected in a theory with an unscreened attractive electro- 
magnetic force between fermion and antifermion. 

This paper will also pursue the more theoretical goal of relating features of 
the theory's spectrum to universal quantities characterizing its renormaliza- 
tion group fixed point. The reader should recall that the continuum model 
must be parametrized by two coupling constants in order for an analytic 
renormalization group to exist |Q, [Q, 0. The electrodynamic coupling 
a must be supplemented by an induced four Fermi coupling G = GqA^/tt^ 
where Go is the dimensionful four Fermi coupling and A is an ultra-violet 
momentum-space cutoff. Analytical studies of the theory's Schwinger-Dyson 
equation reveal that there is a fixed line, in the renormalization group sense. 



extending from a = to a^ = 7r/3 at G = (1 + y'l - a/a^y 0, 0- The 
theory has been particularly well studied at the endpoints of this line. At 
a = 0, G = Gc = 4: the theory reduces to the Nambu Jona-Lasinio (NJL) 



model [^. The cutoff theory brealcs chiral symmetry for all G > Gc and 
when the cutoff is removed the theory becomes a free field. At all other 
points along the fixed line there is an unscreened massless photon which 
guarantees that the critical theory is interacting, with a rich mass spectrum. 
At the other end of the fixed line ( G = 1, a = 7r/3), there is the Miransky- 
Bardeen-Leung-Love point where "collapse of the wavefunction" occurs ^, 
[^]. Critical indices along the line of chiral transitions vary continuously as 
a varies from zero to vr/S 



ly 



l/2y^l - a/a, 

(3mag = (2 - ^1 - a/a,)/2sjl - a/a, (1) 

5 = (2 + V'l - a/a,)/{2 - ^1 - a/a,) 
7 = 1 

with familiar mean field values at a = and Miransky's essential singularities 

at a = ttc [0], H- 

A major objective of any numerical lattice approach to this model would 
be a verification of the continuum model's phase diagram in the (a, G) plane. 
Lattice studies which tune the Action and move in a controlled fashion along 
the fixed line have not yet been done. Instead, non-compact quenched QED 
has been simulated by generating photon field configurations using a fast 
Fourier transform and then calculating fermion propagators in that back- 



ground using staggered lattice fermion fields |T2|. It has been argued that 
the discreteness of the Lattice Action necessarily generates some four Fermi 
couplings, so when the lattice fine structure is increased and a chiral tran- 
sition is found, one is actually simulating a continuum Action with both 
photon exchange and four-Fermi forces. Large lattice, high statistics studies 
have measured the critical lattice coupling /3 = 1/e^ = .257(2) and critical 
indices 5 = 2.2(1), /imag = .78(8). These values match the indices b and /imag 
in Eq.(l) if {a,G) = (.44ac,3.06) [jl3|. These rather precise measurements 
use the framework of the chiral Equation of State (EOS) as has been dis- 
cussed in detail before and will be reviewed below. Briefly, this approach 
requires that we simulate the model within its scaling window but not nec- 
essarily directly at its critical point. Simulation data near but not at (3, and 
near but not at zero fermion mass m should satisfy a chiral EOS and lie on 
a universal curve if "reduced" quantities are plotted. This approach, which 



is borrowed from ancient studies of ferromagnetic transitions, has proved to 
be very successful. 

The model's scaling laws and critical indices can also be found from its 
spectrum of states. The development and exploitation of relations between 
spectroscopy, scaling laws and critical indices is a major theme of this work 
and new analytic and numerical results will be presented. Our most success- 
ful result will be a scaling law between the chiral condensate < ipif) > and 
any dynamically generated mass Ma, 

< ijij >= CaM^^'^^''' (2) 

We shall argue that Eq.(2) should hold anywhere in the theory's scaling 
region, i.e. for couplings near the critical point and symmetry breaking fields, 
such as the bare fermion mass, which are sufficiently weak. The theoretical 
assumptions underlying Eq.(2) consist of (1.) validity of the chiral Equation 
of State and (2.) hyperscaling, the idea that the theory's critical singularities 
are due to the divergence of a single length scale, the correlation length. 
Using our accurate < ^ipip > data on 16^ and 24^ lattices and results for Mp 
on 32 X 16^ lattices we shall find that Eq.(2) is well satisfied with the critical 
indices 

Pmag/l^ = 1.25(1) (3) 

which is in excellent agreement with Eq.(l) which also gives 1.25 at the point 
{a, G) = {A4ac, 3.06) picked out by the lattice Action. In mean field theory, 
such as the NJL model, this ratio of critical indices is unity. Thus, the lattice 
simulation is correctly accounting for the effects of an unscreened vector force 
which renders the theory non-trivial in the continuum limit. 

Other analyses discussed below will give further evidence that the critical 
index 7, which controls the susceptibility divergence, is unity, and u, which 
controls the correlation length, is 0.68(2). 

This article is organized into several additional sections. In Sec. 2 we 
review the relation between the chiral Equation of State, hyperscaling and 
critical indices. The approach and results in ref. play a crucial role here. 
In Sec. 3 we present our numerical results which are then used in Sec 4-6 
to extract the physics discussed briefly above. Finally, in Sec. 7 we suggest 
future research. 



2 Scaling and critical indices 

In this section we shall summarize the main theoretical arguments which 
proved useful in our analysis of our quenched data. A detailed exposition 
can be found in refs. [^, [|14|. The basic philosophy of refs. 0, [1^ is the 



use of the chiral Equation of State and hyperscaling, the assumption that 
a theory's critical singularities are all due to one divergent length scale, to 
obtain scaling laws for low energy quantities which have an immediate physi- 
cal significance. This approach will yield new ways to extract critical indices 
for low energy properties of the theory, and will shed interesting insights on 
these two seemingly different features of the theory. 

We begin with a theory's Equation of State (EOS) written in a standard 



form [|15[. Let t denote the deviation from the critical coupling, and m 
the symmetry breaking field; in our application t = (3 — j3c where j3 = 1/e^ 
( e = electrodynamic coupling) and m is the bare fermion mass which ex- 
plicitly breaks the continuous chiral symmetry of massless QED. The EOS 
records the order parameter's response to a change in coupling and symmetry 
breaking field, 

m =< '4)11) >^ f {t / < l/jlp >V/3ma9^ ^4^ 

Eq. (4) contains the critical indices 6 and Pmag in standard statistical me- 
chanics notation. Eq. (4) should hold everywhere in the theory's scaling 
region, small t and m. It is frequently convenient to rewrite Eq.(4) in terms 
of "reduced" variables and to invert it. Let yhe a reduced symmetry-breaking 
field and 2; be a reduced order parameter 

y = m/t 

(5) 

where we have used the conventional definition A = 5(3mag- If "we divide Eq. 
(4) by t^ we have 

y = z^f{z'^/P^-<>) (6) 

In other words, the EOS implies that the reduced symmetry-breaking 
field is just a function of the reduced order parameter. Upon inverting this 
dependence, 

< ^V' >= CagF{y) (7) 



which is an alternative and famihar form of the EOS. Eqs. (4) and (7) have 
been used extensively to extract the critical indices and the critical coupling 
in quenched lattice QED simulations |]r3|. They will be used again with 
impressive success in Sec. 4 below. 

Our next goal is the determination of "Equations of State" for other low 
energy features of the model, such as its mass spectrum. In order to do that 
we assume the correlation length scaling idea, which underlies hyperscaling, 
familiar in classical statistical mechanics. We assume that there is a single 
macroscopic length scale ^ whose divergence is responsible for the theory's 
critical behaviour. ^ sets the scale for low energy physics (low on the scale 
of the cutoff, the reciprocal of the lattice spacing in our case) and one can 
do dimensional analysis using ^ as a unit of length. In particular, if the field 

has scale dimension d^ and we consider the scaling region of the theory 
where ^ is much larger than the lattice spacing, then 

< >oc r"^^ (8) 

Futhermore, suppose that the theory has power-behaved critical singularities, 
so ^ diverges as t~^ in the scaling region. Any dynamically generated mass 
scale Ma in the model must, by the assumption of correlation length scaling, 
behave as ^~^. So, Eqs. (7) and (8) give, 

< >(x Mi^ oc t^"^ (X tl,,F{y) (9) 

which implies that d^ = Pmag/^ and, 

< >= C^Mf™"/'^ (10) 

everywhere in the scaling region. The constants Ca are not controlled by 
this scaling argument and are non-universal. We shall find that numerical 
data for < iJjiIj > and Mp satisfy the scaling law Eq.(lO) beautifully and 
yield a relatively precise estimate for Pmag/^ which agrees with continuum 
calculations based on the Schwinger- Dyson equation |0], at the point 
(.44ac, 3.06) picked out previously by the lattice results. 

Since the anomalous dimension rj of the field is defined as d^ = d/2 — 

1 +ri/2, Eq.(lO) is directly sensitive to rj itself. The anomalous dimension r) 
vanishes identically in mean field theory. Therefore, Eq.(lO) will prove to be 
particularly informative in this and, hopefully, other theories studied in the 
future. 
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We shall also find it useful to write Eq.(9) in the form 

Ma = f'GM (11) 

and consider ratios of dynamically generated masses, 

R,, = MJM, = Gab{y) (12) 

which are functions of just one scaling variable y. The validity of this scaling 
law will also be demonstrated below. In fact, in the following Sections of 
this article the three equations Eqs. (7), (10) and (11) will be tested and 
used to extract critical indices. The vahdity of hyperscaling will be tested 
by plotting two mass ratios against one another. The universality of these 
plots will provide an independent check of the scaling hypothesis which was 
first seen in EOS plots of < ipif) >. The limiting (m — ^ 0, t — *■ 0) values of 
particular mass ratios are universal numbers and we shall attempt to estimate 
them below. 

In ref fl^ the properties of the mass ratios are discussed in great detail. 



and several analytic examples are given. It is observed there that in a wide 
class of theories in which mesons are made out of fundamental, pointlike 
constituents, non-triviality and compositeness go hand-in-hand. In particu- 
lar, mesons have non-zero sizes which guarantee that they interact with one 
another. Only in the mean field limit where the size of mesons shrinks to 
zero and compositeness is lost, do the theories become trivial. In that case 
the mass ratio 2Mf/M„ is exactly unity. We shall find in quenched QED 
that Mcr < 2Mf in the chiral limit, which is indicative of compositeness and 
non triviality. 

We direct the reader to ref. [Q , [^ for a more thorough discussion of the 
connection between hyperscaling, renormalizability and non-triviality. We 
have simply extracted a few portions of a wide class of relations for our use 
here and have emphasized only those which shed light on the non-triviality 
of the light quark dynamics in quenched lattice QED. 

3 Numerical Results 

We describe here the numerical simulations and the data analysis. The 
reader uninterested in these details can skip this section and go to the sub- 
sequent sections devoted to the discussion of our results. 
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3.1 The Simulations 

We have collected 100 , 30 and 248 gauge field configurations on the IG'*, 
24^ and 16^ x 32 lattices, respectively , by using the algorithm introduced in 
[|T^ . The algorithm begins in momentum space and produces the appropriate 
Gaussian distribution of photons. Using a Fast Fourier Transform it then 
generates a set of dimensionless gauge fields 6^^ in coordinate space. The 
coupling of the gauge field to the fermion is implemented by rescaling the 
gauge fields as 6* — >• B\f^. 

This fast algorithm eliminates the correlations between subsequent con- 
figurations (the only source of correlation comes from the random number 
generator) and greatly reduces the correlations between the fluctuations at 
different spatial sites on the lattice. Since this algorithm produces indepen- 
dent configurations, our statistical sample is quite large by present simulation 
standards. 

3.2 The Chiral Condensate 

The configurations produced on the symmetric lattices have been used for 
the computation of < ■?/'■?/' >. We have inverted the staggered Dirac operator 
for each value of the fermion mass using a second order conjugate gradient 
routine. We used a noisy source for this inversion [|l^. On the 16*^ lattice 
at each /? value (spaced by A/3 = 0.002 and ranging from /3 = 0.247 to 
/3 = 0.257) we have accurate < ■?/''?/' > values at five fermion masses m equi- 
spaced from .001 to .005. We have seen in past studies that accurate data 
over this mass range can yield unambiguous, quantitative results for /3c, (5 
etc. We also recall from past studies that the finite size effects, as inferred 
from 24^ and 32^ simulations, are very small in < ^Z)^ > over this range of m. 
This is further demonstrated by our new data on the 24^ lattice (presented 
in Table 0), where we explored the same set of masses m and a very large 
set of /3's. 

3.3 The Analysis of the Spectrum Data 

The spectroscopy data are from the 16^ x 32 lattice where a wall source 
was used for the inversion of the Dirac operator after the appropriate gauge 
fixing. 



<1p1p > 




0.001 0.002 0.003 


0.004 0.005 


0.257 


0.0392(9) 0.0541(7) 0.0656(7) 


0.0752(7) 0.0833(6) 


0.255 


0.0485(11) 0.0622(8) 0.0736(8) 


0.0828(7) 0.0907(7) 


0.253 


0.0583(12) 0.0716(9) 0.0822(9) 


0.0909(7) 0.0984(7) 


0.251 


0.0694(14) 0.0817(11) 0.0913(10) 


0.0994(8) 0.1065(8) 


0.249 


0.0814(16) 0.0920(13) 0.1006(12) 


0.1081(9) 0.1149(8) 



Table 1: < ipip > data at various fermion masses and couplings (3 on the 16^ 
lattice 

We worked at five (3 values, ranging from 0.245 to 0.265. For each jS we 
have data for five masses, spaced by 0.001, and ranging from 0.002 to 0.006. 
We have thus two sets of data in the weak coupling region (at (3 = 0.260 
and 0.265 ) where we expect to observe clear signatures of chiral symmetry 
restoration. 

For hadron operators we use the standard local form. The propaga- 
tors in the PSeudoscalar, SCalar, VecTor, Pseudo Vector channels are then 



parametrized (in the region of large time) as follows [18 



G(r) 



a[exp{-MT) + exp{-M{Nr - r)] + 
{-iya[exp{-MT) + exp{-M{Nr - r)] 



Turning to the fermion, its parametrization reads: 

G{t) = a[exp{-MT) - {-iyexp{~M{Nr - r)] 



(13) 



(14) 



By fitting to the above forms, we obtained good estimates of the fun- 
damental particles (i.e. the lower masses in the direct channel for the PS , 
VT sectors, and in the oscillating channel for the SC, PV sectors). To keep 
the (possible) contaminations from excited states under control we vary the 
extrema of our fits: typically we fit for ti < t < T/2 (2 < ti < 10), the 
number of degrees of freedom ranging from 6 to 12. The results level-off for 
ti ~ 7 — 8 for the pion and the fermion (for which we performed 1-particle 
fits), and for ti ^ 2 — 3 for the other particles (2-particles fits). In the first 
case (pion and fermion) the errors on the mass estimates were obtained from 
a jack-knife analysis performed by decimating one configuration at a time 
(remember that our configurations are independent). For the other mesons. 



9 



< VV" > 




0.001 




0.002 




0.003 




0.004 




0.005 


0.260 


0.0296 


; 14 


) 0.0436 


; 12 


) 0.0546 


; 11 


) 0.0639 


; 10 


) 0.0721( 9) 


0.259 


0.0324 


! 15 


) 0.0466 


; 13 


) 0.0577 


; 12 


) 0.0670 


; 11 


) 0.0751( 10) 


0.258 


0.0352 


; 16 


) 0.0498 


; 14 


) 0.0609 


; 12 


) 0.0702 


; 11 


) 0.0783( 10) 


0.257 


0.0383 


; 18 


) 0.0532 


; 15 


) 0.0643 


; 13 


) 0.0735 


; 12 


) 0.0816( 10) 


0.256 


0.0419 


;2o 


) 0.0569 


; 16 


) 0.0678 


; 14 


) 0.0770 


; 12 


) 0.0850( 11) 


0.255 


0.0460 


; 23 


) 0.0608 


; 17 


) 0.0715 


! 14 


) 0.0806 


; 12 


) 0.0885( 11) 


0.254 


0.0508 


[ 25 


) 0.0648 


; 18 


) 0.0753 


[ 15 


) 0.0843 


[ 13 


) 0.0922( 11) 


0.253 


0.0559 


; 28 


) 0.0690 


[ 19 


) 0.0793 


[ 15 


) 0.0881 


; 13 


) 0.0959( 12) 


0.252 


0.0611 


; 30 


) 0.0732 


; 20 


) 0.0834 


; 16 


) 0.0921 


[ 14 


) 0.0997( 12) 


0.251 


0.0662 


; 31 


) 0.0776 


; 20 


) 0.0876 


; 16 


) 0.0961 


[ 14 


) 0.1037( 12) 


0.250 


0.0710 


;3i 


) 0.0821 


;2o 


) 0.0919 


; 17 


) 0.1003 


; 14 


) 0.1077( 13) 


0.249 


0.0759 


; 31 


) 0.0869 


; 21 


) 0.0964 


; 17 


) 0.1046 


[ 15 


) 0.1118( 13) 


0.248 


0.0812 


;3i 


) 0.0918 


; 21 


) 0.1010 


; 17 


) 0.1090 


; 15 


) 0.1161( 14) 


0.247 


0.0869 


;3i 


) 0.0970 


; 22 


) 0.1057 


; 18 


) 0.1135 


; 16 


) 0.1204( 14) 


0.246 


0.0930 


;3i 


) 0.1022 


;23 


) 0.1106 


; 19 


) 0.1180 


; 16 


) 0.1247( 14) 


0.245 


0.0990 


; 31 


) 0.1075 


; 24 


) 0.1155 


; 20 


) 0.1227 


; 17 


) 0.1292( 15) 


0.244 


0.1048 


; 33 


) 0.1128 


[ 25 


) 0.1204 


; 21 


) 0.1274 


[ 18 


) 0.1337( 16) 


0.243 


0.1104 


; 34 


) 0.1182 


; 26 


) 0.1254 


[ 21 


) 0.1322 


[ 18 


) 0.1383( 16) 


0.242 


0.1159 


; 35 


) 0.1236 


; 27 


) 0.1305 


; 22 


) 0.1370 


[ 19 


) 0.1430( 17) 


0.241 


0.1215 


; 35 


) 0.1291 


; 27 


) 0.1358 


; 23 


) 0.1420 


; 20 


) 0.1478( 18) 


0.240 


0.1279 


;35 


) 0.1347 


;29 


) 0.1410 


; 24 


) 0.1470 


; 21 


) 0.1527( 18) 


0.239 


0.1344 


;35 


) 0.1406 


;3o 


) 0.1464 


;25 


) 0.1521 


; 22 


) 0.1575( 19) 


0.238 


0.1411 


;37 


) 0.1466 


;32 


) 0.1519 


;27 


) 0.1573 


;23 


) 0.1625( 20) 


0.237 


0.1478 


;4o 


) 0.1526 


;34 


) 0.1574 


;28 


) 0.1625 


; 24 


) 0.1675( 21) 


0.236 


0.1542 


;43 


) 0.1586 


;36 


) 0.1630 


;29 


) 0.1678 


;25 


) 0.1725( 22) 


0.235 


0.1602 


;47 


) 0.1645 


;38 


) 0.1686 


;3i 


) 0.1731 


;26 


) 0.1776( 22) 



Table 2: < i/jip > data at various fermion masses and couplings (3 on the 24"^ 
lattice 
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we quote the weighted average of the results obtained with different ti. We 
also performed effective 2-niass analyses by exactly solving the equations for 
masses and amplitudes on four subsequent time slices; the results are in full 
agreement with the ones we quoted, but the errors are larger. In any case the 
agreement of the results obtained by local and global fits is a nice consistency 
check on our analyses. 

To get a good fit the inclusion of secondary particles was necessary in 
the VT, SC,and PV sectors, as already said, but their amplitudes proved 
to be very small, and often compatible with zero considering the statistical 
errors. In many cases we noticed that the inclusion of a secondary state 
greatly improved the quality of the fit (i.e. reducing the x^) while leaving 
completely unaffected the results for the fundamental states. In some cases 
the results of the fits suggest the existence of a very light particle in the 
secondary channel. We confirmed that the numerical significance of very 
light states was less in the present data as compared to past simulations 
on smaller lattices. We believe that any evidence for a very light state is 
unreliable and is, in fact, indicative of rather large finite size effects. 

Figs. 1-3 provide an overview of the results of the fits: we note from them 
that the relative amplitude of the oscillating channel to the non-oscillating 
one decreases while going from weak to strong coupling, and from low to 
high masses. Figs. 1 shows a sample of vr propagators for various sets of 
parameters, with the best fits superimposed. Fig. 2 contains several typical 
fits for the vector channel. We plot there the fit obtained with t = 5 as a 
starting point. Note that the fit predicts the data at small t (i.e. t < 5.) 
with good accuracy. Some results for the the sigma are shown in Fig. 3 
: its behaviour is satisfactory in the strong coupling region. In the weak 
coupling region the situation is somewhat less controlled, mostly because of 
large errors in the large t region, and we feel that the results for a at weak 
coupling should be considered with some extra care. 

Summarizing, the weak coupling, low mass region turns out to be rather 
difficult to treat, while the results for the fundamental particles in the critical 
and strongly coupled region are fully satisfactory. Tables 1^-0 collect our 
results. 
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TT 




0.265 0.260 0.255 0.250 0.245 


0.002 


0.0909( 17) 0.1090( 18) 0.1273( 18) 0.1419( 18) 0.1479( 16) 


0.003 


0.1189( 13) 0.1402( 16) 0.1589( 16) 0.1745( 18) 0.1794( 13) 


0.004 


0.1458( 13) 0.1668( 13) 0.1871( 16) 0.2012( 16) 0.2070( 15) 


0.005 


0.1685( 13) 0.1917( 15) 0.2109( 15) 0.2243( 15) 0.2298( 13) 


0.006 


0.1911( 13) 0.2133( 13) 0.2322( 13) 0.2450( 16) 0.2500( 12) 



Table 3: Results for the pion mass from one particle fit 



p 




0.265 0.260 0.255 0.250 0.245 


0.002 


0.1653( 29) 0.2044( 42) 0.2520( 51) 0.3231( 83) 0.3934(118) 


0.003 


0.1916( 26) 0.2310( 29) 0.2815( 37) 0.3472( 70) 0.4164( 86) 


0.004 


0.2163( 23) 0.2579( 27) 0.3089( 41) 0.3699( 52) 0.4355( 71) 


0.005 


0.2385( 23) 0.2797( 22) 0.3322( 39) 0.3919( 40) 0.4548( 62) 


0.006 


0.2598( 19) 0.3026( 22) 0.3535( 35) 0.4115( 41) 0.4734( 59) 



Table 4: Rho masses from a two particle fit 



a 




0.265 0.260 0.255 0.250 0.245 


0.002 


0.1905(220) 0.2223(403) 0.2361(259) 0.2922(282) 0.3418(691) 


0.003 


0.2182(132) 0.2338(111) 0.2681(427) 0.3194(281) 0.3690(536) 


0.004 


0.2287( 82) 0.2625(427) 0.2959(200) 0.3480(410) 0.4028(276) 


0.005 


0.2469( 64) 0.2796(111) 0.3207(184) 0.3801(529) 0.4271(361) 


0.006 


0.2673( 80) 0.3032( 82) 0.3447(137) 0.3950(225) 0.4531(170) 



Table 5: Results for the sigma mass from a two particle fit 



al 




0.265 0.260 0.255 0.250 0.245 


0.002 


0.1829( 38) 0.2151( 44) 0.2692( 89) 0.3434(425) 0.3716(135) 


0.003 


0.2057( 31) 0.2441( 50) 0.2850( 40) 0.3638(216) 0.4051(103) 


0.004 


0.2271( 26) 0.2601( 43) 0.3088( 38) 0.3767( 70) 0.4307( 83) 


0.005 


0.2463( 26) 0.2837( 34) 0.3315( 37) 0.3950( 74) 0.4544( 82) 


0.006 


0.2682( 44) 0.3043( 30) 0.3528( 38) 0.4127(138) 0.4741(119) 



Table 6: al masses from a two particle fit 
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fermion 




0.265 0.260 0.255 0.250 0.245 


0.002 


0.0951( 50) 0.1116( 52) 0.1293( 92) 0.1481(151) 0.1782(231) 


0.003 


0.1071( 40) 0.1244( 44) 0.1375( 96) 0.1623(134) 0.1905(214) 


0.004 


0.1180( 37) 0.1343( 61) 0.1497( 89) 0.1742(122) 0.2022(179) 


0.005 


0.1283( 34) 0.1407( 68) 0.1607( 85) 0.1857(116) 0.2136(166) 


0.006 


0.1383( 33) 0.1506( 66) 0.1714( 83) 0.1963(111) 0.2244(157) 



Table 7: fermion masses, from one particle fit 

3.4 The Chiral Extrapolation 

We conclude this section by briefly describing the extrapolation procedure 
which leads to the values quoted in Table || and plotted in Fig.4. These tech- 
nical points are important since they underlie the evidence for the Goldstone 
behavior of the pion we will discuss in more detail in the next Section. 

A crucial ingredient for the chiral extrapolation is the knowledge of the 
functional dependence of the hadron masses on the bare quark mass. This is 
not known a priori, and could depend on f3, possibly being sensitive to the 
anomalous dimension developed while approaching the critical point. In an 
attempt to use as little prejudice as possible, we decided to fit our data to 
several different functions of m. We considered quadratic and cubic hadron 
masses as first and second order polynomials in m. We also tried a Fade' 
P[l,l] expression M = [a + bm)/{l + cm). To critically evaluate the approach 
to zero of the extrapolated values as the coupling approaches the critical point 
(or its consistency with zero in the case of Pion), we tried fits both with and 
without a constant term. 

All the fits clearly favour the standard form M^ = MQ+Am + Bm'^ for all 
the mesons. Mq is zero for the pion, and non-zero for the other particles in 
the strong coupling region. The extrapolated data are not quite zero in the 
weak coupling region as well because, we beheve, of finite size effects which 
push the masses up. Non-zero extrapolated values are expected (and found) 
in the weak coupling region for the same reason. It is interesting, and gives 
us more confidence in the results, that the Mq values obtained with different 
parametrizations are quite similar to one other. 

In Fig.4 we present an overview of the results: (fermion, a , p, al) {m = 0) 
in the strong coupling region are plotted vs. /?. In this plot we discarded 
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mesons/fermion (m = 0) 




p al (T f 


0.245 


0.3466( 14) 0.2867(101) 0.2623( 58) 0.1513( 6) 


0.250 


0.2664( 24) 0.3207( 26) 0.2190( 39) 0.1150( 25) 


0.255 


0.1721(105) 0.2274( 50) 0.1547(111) 0.1119( 25) 


0.260 


0.1367( 64) 0.1611( 36) 0.1895(125) 0.0785( 58) 


0.265 


0.0965( 84) 0.1337( 28) 0.1694( 15) 0.0671( 20) 



Table 8: Extrapolated values according to M^ 



Mn + am + 6m^ 



the data for the al at the lowest beta value. It appeared that M ~ 0.4 is an 
upper limit to the masses one can meaningfully study at /3 = 0.245 because 
of finite spacing effects. The finite size effects show up clearly in the results 
for the fermion mass near (3c'- the fermion masses "flatten" near the critical 
point to a value around 0.1 . Indeed 1/0.1 is close to half the time extent of 
the lattice, so it is unlikely that we can reliably calculate lower masses. 

All the extrapolated masses are well split in the strong coupling region 
and approach zero (with the caveats discussed above) as /? — > Pc- 

4 The Equation of State and the Critical In- 
dices 7, S 

Our previous work on chiral symmetry breaking in quenched QED4 has 
shown that its critical behaviour is not described by mean field theory [|13 



In fact, our measurements have all been consistent with the quantitative 
predictions of the Schwinger- Dyson equation defined with a momentum cutoff 
A and the hyperscaling relations between critical indices 0, [§!• Our purpose 
here is to confirm those results using new, more accurate data and to extract 
the critical index 7, which controls the susceptibility divergence in the critical 
region. A practical way to find 7 proceeds through the chiral equation of 
state, 

<ipi;> lm>l^ = f{A(3/ < ipi^ >i//3-.9) (15) 

where 6 and Pmag are the usual critical indices discussed in Sec. 2 above, 
A/? = P — Pc and / is a universal function. Our past low mass < ipip > data 
has given j3c = 0.257(1), 6 = 2.2(1) and Pmag = .78(10). These results were 
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also obtained by combining Lanczos data for the spectral density of < i/jiIj > 
with conjugate gradient data for its bare fermion mass dependence. Here 
we used the conjugate gradient method to calculate < ipip > and FFT's to 
generate independent background photon configurations |[16[| . Our < ipip > 



data is presented for 100 independent gauge field configurations on a 16^ 
lattice in Table |I|. 

The data in Table ^ picks out Pc = 0.257(1) very clearly. Recall that for 
P = Pc the order parameter < ipip > should scale as 

< iji/j >= Am^'\p = (3c) (16) 

We test for a simple power law by plotting log < ■?/'-?/^ > vs. logm for (3 = 0.257 
and 0.255 in Fig. 5. Note that the f3 = 0.255 data indicates that < ipip > 
at m = 0.001 is too large to admit a good power law fit for m ranging from 
0.005 to 0.001. However, the power-law fit for the /3 = 0.257 data is excellent 
and the slope of the line in Fig. 5 yields 6 = 2.15(5) in agreement with our 
earlier work. Inserting these values for f3c and 6 into the chiral equation 
of state (1) we can ask whether there is a particular Pmag which yields a 
universal curve /, i.e. can all the data of Table ^ at various (3 and m values 
be plotted on a single scaling curve? Recalling the definition 7 = Pmagi^ — 1) 
then gives us the susceptibility index of interest. In Fig. 6 we show the chiral 
equation of state for /3c = .257, 6 = 2.20 and Pmag = -833 The success of 
scaling hypothesis is very impressive. This value of Pmag yields 7 precisely 
equal to 1.00. This result is consistent with, but more accurate than, our 
past analyses and is in agreement with the Schwinger- Dyson prediction of 
7 = 1.00. We discussed the theoretical significance of this result in Sec. 2 
above. 

The almost perfect scaling seen in Fig. 6 deteriorates rapidly as 7 devi- 
ates from unity. The dispersion in / exceeds the statistical errors (several 
representative error bars are shown in Fig. 6a) if 7 exceeds 1.05 or if 7 falls 
below .975, so we have determined 

7 = 1.00 + 0.05-0.03 (17) 

We can exploit the scaling form of EOS also by writing 

m/ < tPij >^= f{{(3c - 13)/ < # >^/^"''"0 (18) 
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On the basis of the very detailed analysis of ref. [0 we can predict that the 
universal function / is straight line. In fact we know that the data for the 
chiral condensate satisfy 

m = a(/?c - P) < ^p'^p > +b < ipi; >^ (19) 

which confirms the theoretical expectation 7 = 1. The previous equation can 
be restated like this 

m/ < ipip >^= a{(3 - /?c) < i^i! >^~^ +b (20) 

We can then substitute into the EOS 

a{(3 - 13c) < iJifj >^-^ +b = f{{f3c - (3)/ < ipip >^/^'-'^^) (21) 

So f must be a linear function with l/Pmag = ^ — S, i-e. 7 = 1.0. It is 
interesting to check that / actually is a straight line by looking at the data 
on the 24"^ lattice, shown in Table |^, which explores a wide set of /3's. The 
linear behaviour of Fig. 6b is very clear, and the conclusion 7 = 1.0 receives 
further support. 

5 Scaling and the Critical Indices 77 and u 

As discussed in Sec. 2, we can obtain estimates of u and rj by exploiting 
the relation 

< ^^ >= C'M^°'/2-l+r,/2) _ Q^j{f3^a,M (22) 

which should hold everywhere in the scaling region. Moreover, we can test the 
consistency of u by plotting Mp/t" vs m/t^'""^'^ as suggested by the equation 
of state for the masses 

Mp = r^(m/t'^'""«'^) (23) 

We will see that the estimates of u obtained by using the EOS are by far 
more convincing that the one based on the direct fit at m = 

Mp = A{P - P,y (24) 

which is however consistent with the previous ones. 
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5.1 < i/ji/j > vs Mp in the Scaling Region 

Consider the scaling law, 

< 2pij >oc M^^™'"^/^ (25) 

reviewed in Sec. 2 above. Since mean field theory predicts that Pmag = ^i 
deviations from linear dependence of < ■?/'■?/' > on Mp provide clear evidence 
of a non-trivial critical point. The Schwinger- Dyson prediction for the con- 
tinuum model is (imaglv = 1.25 if 6 = 2.2, as determined from our Equation 
of State fit in Sec. 3. 

We took < ipilj > data from the 24^ run. They overlap with the spec- 
troscopy data at /9 = ( 0.260, 0.255, 0.250 , 0.245), m = (0.002, 0.003, 0.004, 
0.005). Note that f3 = 0.260 lies in the weak coupling region. However since 
the transition is rounded, the data obtained at this [3 value can be (tenta- 
tively) added to the data in the strong coupling region . We analyzed data 
sets both with and without the (3 = 0.260 points and confirmed that our fits 
did not change much. 

Fig. 7 shows a log-log plot of < ipip > versus log (Mp), All the points lie 
on the same universal plot (The only small deviation is at /? = 0.260, at the 
smallest mass value.) We stress that this is an highly non-trivial test of Eq. 
(25), considering that the data we used are from completely independent 
simulations. The straight line is our best fit with all the data included : 



1.25 

■ _ _ 

We fit both Mp vs < ipip > and < ipip > vs Mp, with and with- 



< ^^ >= Mp 



out the points at /9 = 0.260. The four fits give, respectively i^maglv = 
1.265, 1.275, 1.243, 1.249. Assuming that Pmag = 0.86(3) we get for u: u = 
0.680(25), 0.674(25), 0.691(25), 0.688(25). These values are different from 
the mean field prediction, u = 1/2, but it is not easy to estimate the errors 
associated with this fitting procedure. Rather than attempting this, we pre- 
fer to show the inconsistency of our results with the mean field prediction in 
a more direct fashion. In Fig. 8 we plot < ipip >^ vs Mp for different x values. 
The solid line corresponds to our best fit, the dashed ones to the mean field 
result X = 1, the dotted line to Pmag = 0.86 and u = 0.5. The sensitivity of 
the quality of the fits to the choice of the exponent x is clear, and thus the 
incompatibility of our results with the mean field prediction is established. 

We see that the anomalous dimension rj which follows is ~ 0.5. This is 
a remarkably large value. Nonetheless, it is consistent with the hyperscaling 
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relation between 6 and 77: 

(6-r7)/(2 + r/) = 5 (26) 

with S = 2.2 as determined in Sec. 3 above. 

In conclusion, our data rule out mean field behaviour. They predict 
1/ = 0.68(3) if /5^,g = 0.86(3). 

5.2 The Equation of State for the Masses and Corre- 
lation Length ScaUng 

As stated in Sec. 2, we can check our estimate of u by using the data of 
the masses alone. This is done in Fig. 9 where we plot Mp/{Pc — PY vs 
m/{(3c - /?)^-''«^ for u = 0.67, p^nag = 0.833, 5 = 2.20. We use only points at 
strong coupling. The scaling hypothesis works very well. Again, we checked 
that the scaling behavior deteriorates if we select mean field exponents. The 
equation of state for the masses can be studied for m = 0. In this case 
we get Mpijn = 0) = g{S^)t" and a direct measure of z/ can be obtained 
by fitting Mp{m = 0) as a power of /5 — [3c. Unfortunately the quality 
of our data does not allow a careful estimate of v in this way because of 
uncertainties in the extrapolation to m = 0. The point at /? = 0.255 can 
be affected by finite volume effects, and the two remaining points alone are 
not sufficient to allow a meaningful, controlled estimate of v. It is necessary 
to work very close to the transition, because previous studies have shown 
that in the strong coupling region mean field theory functional dependences 
are expected. With these caveats it is still interesting to get an independent 
estimate of v even if it is crude. Our best fit (which is shown in Fig. 10) 
gives Mp{m = 0) = A{f] — 0.260)^'^"^^ which is nicely consistent with the more 
sophisticated analysis done in Sec. 5.1 above. It also interesting to note that 
the ratio of the extrapolated masses for the p and the fermion is constant 
within statistical errors: we obtained Mf/Mp = 0.434(5) at /5 = 0.245, and 
Mf/Mp = 0.426(2), at /? = 0.250, in reasonable mutual agreement. This 
agreement again supports scaling behaviour. 

5.3 Critical indices, summary 

We present in Table ^ an overview of the critical exponents and the 
hyperscaling relations. Let us recall that the estimates for Pmag, S, rj, v 



Critical indices 


Index 


Result from the simulation Result from HS 


MF 


Pmag 


0.86(3) 0.86(6) 


0.5 


1 = (3mag{S -I) 


1.00(5) 


1.0 


6 


2.2(1) 2.16 


3.0 


V 


0.5(1) 0.5 


0.0 


V 


0.675 0.68(3) 


0.5 


-AU + 2l3mag + 7 


0.1(1) 0.0 


0.0 


{2-r^)vh 


1.1(1) 1.0 


1.0 


{2v-^)/{vr,) 


1.1(1) 1.0 


1.0 


{6-rj)/{2 + r^)-6 


0.13(20) 0.0 


0.0 



Table 9: Critical indices and relations among them: in the first column results 
from the simulation, in the second column the HS(hyperscaling) prediction 
assuming the other indices as input for the single index entries, in the third 
column the mean field prediction (d=4) 

come directly from independent numerical analysis, and 7 has been computed 
according to 7 = Pmagi^ ~ !)• So, in the first column we quote the results of 
the numerical analysis, and in the second column the results obtained from 
the hyperscaling relations assuming the other exponents as input. (When we 
do not have a safe estimate of the errors we simply quote the central values.) 
The third column shows the mean field results (we recall that 7 = 1.0 is also 
the theoretical prediction). We summarize in the last four rows of Table ^ 
the hyperscaling relations we used. 

We believe that these results support a picture of non-trivial critical be- 
haviour consistent with hyperscaling, with a large anomalous dimension 77. 



6 Spectroscopy: Results 

We devote this Section to the discussion of the Spectroscopy results 
presented in Tables. |]- 0. 
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6.1 The Goldstone Character of the Pion 

Our first task is to check that the pion of quenched QED has all the at- 
tributes of a Nambu- Goldstone particle associated with spontaneous break- 
down of a continuous chiral symmetry. These properties include 1. the pion 
should be massless in the strong coupling phase where chiral symmetry is 
broken, 2. the square of the pion mass should be proportional to the bare 
mass of the fermion if that mass is sufficiently small, and 3. PCAC relations 
should hold for f.„,m,M.„, and < ipip >■ 

Using the data in Table 3 we checked that M^ = Am + Bm^ works well 
(see Section. 2 for more details) : this is the usual parametrization for the 
dependence of the pion mass on the bare quark mass assumed in QCD when 
the quark masses are not really small. We checked that if we added a constant 
term to the fit, its best value was consistent with zero. 

We also looked at the dependence of the pion mass on the quark mass in 
a slightly different way, i.e by fitting M^ = Am^ (Fig- H)- We see that in 
the strong coupling region the slopes of the log-log plots for different /5's are 
quite similar, and roughly consistent with a square root dependence of Myr 
on m (the actual values of the slopes are 0.478, 0.495, 0.547 0.610 0.678 from 
/3 = 0.245 to /3 = 0.265). In conclusion, we observe only slight deviations 
from Mjr ~ J{'m) in the strong coupling region. These deviations can be 
eliminated by letting the exponent vary slightly, or, more conventionally, by 
adding a second order term in the m expansion. In the weak coupling region 
the deviations are somewhat more important, but they can be controlled 
in the same way. In conclusion, the pion has the properties of an ordinary 
Goldstone boson. 

To complete the discussion on PCAC, we show in Table |10] (and in Fig. 12) 
the results for 

fl = 2m < ip^ > /ml 

The data for < i/jiIj > overlap with the spectrum data for four mass values 
and four /3 values. There is sufficient data to see that /^ is nicely linear 
in m (the lines in the plot are to guide the eye), and, most importantly, 
extrapolates to a non-zero value. 
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fl 




0.002 0.003 0.004 0.005 


0.260 
0.255 
0.250 
0.245 


0.0147( 6) 0.0167( 5) 0.0184( 4) 0.0196( 4) 
0.0150( 6) 0.0170( 4) 0.0184( 4) 0.0199( 3) 
0.0163(5) 0.0181(5) 0.0198(4) 0.0214(3) 
0.0197( 6) 0.0215( 5) 0.0229( 4) 0.0245( 4) 



Table 10: /,r as a function of (3 and m 

6.2 Mass Ratios and Level Ordering 

We discuss now the behaviour of some relevant mass ratios: we present 
them in a scale invariant form in Figs. 13-16 , where we plot only the points 
belonging to strong coupling region: the squares are for j3 = 0.255, the 
diamonds for j3 = 0.250, the crosses for 0.245 . Note that since we used 
the same values of quark masses at the different (3 values the physical mass 
region we explore changes with /?: near the critical point the physical masses 
are bigger. 

Let us first discuss the behaviour of p/ fermion. The plot (Fig. 13) shows 
nice scaling (i.e. all the points fall on top of each other). Some deviations 
are observed near the critical point (but still the data points are on the same 
line within errors), which would suggest a p lighter than twice the fermion 
mass in the chiral limit. However, practical difficulties, discussed in Sec. 3 
above, which effect the fits of the fermion propagator at small mass near 
criticality, together with the good agreement of the scaled data at (3 = 0.245 
and /3 = 0.250 , suggest that the leftmost point at /3 = 0.255 be discarded so 
that Mp{m = 0) > 2M/(m = 0). 

The data for a/ fermion (Fig. 14) are particularly interesting for techni- 
color applications [^ and show nice scaling as well (unfortunately the statis- 
tical errors are large). The results in the figure suggest that in the chiral limit 
the sigma mass is less than twice the dynamical mass of the fermion. How- 
ever, the deviations from the NJL result (that the sigma is twice as heavy as 
the fermion) are not very large and are partially masked by statistical errors. 
As discussed in Sec. 2 above, and in [Q a deviation from the NJL result for 
this ratio implies non-trivial behaviour. 

The al and p (Fig. 15) are almost degenerate in the large mass region. 
Some deviations to scaling are observed at /? = 0.245: they can be ascribed 
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both to true scaling violations, and finite spacing effects (note that the mag- 
nitude of finite spacing effects depend on masses in lattice, not in physical 
units). The figure suggests that the degeneracy of the al and the p is resolved 
in the chiral limit, and the al is the heavier state. 

The last plot shows the ratio between the a and p (Fig. 16), and no 
deviations in scaling are observed. 

In conclusion, our data show some evidence for asymptotic scaling ( con- 
tinuum behaviour), and suggest the level ordering < a < 2 x fermion < 
p < al. 

7 Conclusion 

The results presented here on spectroscopy and critical indices of quenched 
lattice QED are in good agreement with analytic calculations based on the 
continuum formulation of the model. The hypothesis of correlation length 
scaling seems well verified by our data which shows universal behaviour wher- 
ever it is expected. Our analysis of the relation < tp^ >oc M^^°-o^'^ was 
particularly precise and gave a value fimag/^ = 1.25(1) in excellent agree- 
ment with our past interpretation of lattice studies of the model's EOS - 
that the lattice Action corresponds to the point (a, G) = {Aiac, 3.06) on the 
renormalization group fixed line of the continuum model. 

Our spectroscopy results illustrate the type of techni-meson masses one 
should expect in a class of models of the top quark. It was particularly inter- 
esting that our data favored the inequality M„ < 2Mf expected of a theory 
with unscreened vector forces rather that the familiar NJL result that the 
techni-sigma lies at the two fermion threshold. Earlier analytical work on 
quenched planar QED had suggested that the critical point would be char- 
acterized by a massless dilaton. This picture is not supported by our data. 
Later analytic work recognized that dilaton current conservation is broken 
by hard operators and the techni-sigma has no reason to be light. "Realistic" 
Technicolor and Extended Technicolor models of the dynamically generated 
mass of the top quark are not left-right symmetric like the lattice model con- 
sidered here. There are proposals in the literature for incorporating chiral 
fermions into lattice regulated theories. If any of these proposals prove the- 
oretically sound, and computationally practical, then it would be interesting 
to generalize the calculations of this paper to other, more realistic models of 
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the dynamics behind the Standard ModeL 

Instead of reiterating the results of Sees. 4-6 above, we close with a 
short discussion of the puzzling features of the quenched model and our 
simulations. Since the underlying photon dynamics is perfectly free, the 
quenched model does not satisfy all the consistency conditions of a bona- 
fide field theory. However, we do not yet clearly understand the model's 
limitations. In fact, we found that its light fermion- meson sector satisfies 
all the scaling laws expected of complete theories. One peculiarity however 
is the observation that the simulations in the weak coupling phase do not 
display all the expected results. For example, since chiral symmetry is not 
broken at weak coupling the pion and the sigma should be degenerate, and 
the p and al also should be degenerate here. No evidence was found for this 
symmetry restoration in our spectroscopy data. This failure could just be 
technical - we observed that the usual fitting procedures were not compelling 
at weak coupling, so other functional forms of the meson propagator's spatial 
dependence should be investigated. Another problem with the data concerns 
the dynamical fermion mass. We found that this mass "flattens" out as the 
bare fermion mass is taken small, so it does not have the same systematic 
behaviour as the meson masses. We suspect that this result is due to the 
fact that the fermion carries an unscreened charge and therefore experiences 
the long range interactions of the free photons. Under these conditions its 
propagator probably suffers from especially large finite size effects which push 
the dynamical mass estimates up when the bare fermion mass becomes too 
small. It would be useful to pursue this point and obtain better control over 
the fermion mass since it plays such a central role in theoretical developments. 

And finally, lattice investigations of quenched QED must be generalized 
to move away from the critical point (q;,G) = (.44q;c, 3.06). We must learn 
how to tune the lattice Action to move along the critical line predicted by 
the analytic calculations. If this cannot be done, then there may be some in- 
gredient missing in the analytic calculations, or some unsuspected limitation 
in the lattice approach. From a physics perspective, it would be particu- 
larly interesting to tune the lattice Action and watch the critical indices vary 
continuously as we move from the free field NJL point, to the super-critical, 
essential singularities of the Miransky point. 
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Figure Captions 

1 a)Pion propagators at /9 = 0.260 as a function of euclidean time for tlie 

five mass values. The solid lines superimposed are the results of the 
fits for 5 < t < 16. Fig. lb is the same as Fig. la, at /3 = 0.250. 

2 a), b) As Figs. 1, for the p. 

3 <i)(j propagators at /9 = 0.250 as a function of euclidean time for m = 0.002. 

The solid line superimposed are the result of the fits for 5 < t < 16. 
Fig. 3b is the same as Fig. 3a, for m = 0.006. 

4 Results for the chiral limit for {fermion, a, p, al) ( circles, diamonds, 

squares, crosses ) (obtained according to M^ = Mq + am + bm'^) as 
a function of f3. 

5 log < tptp > vs logm at P values 0.257 (squares) and /3 = 0.255 (circles). 

Data from table 1. 

6 a) Chiral equation of state for Pc = 0.257, 6 = 2.2 and Pmag = 0.833 on the 

16"^ lattice. {Pc - P)/ < i'i' >-^/'^™''s is plotted versus < ipi/j > /m^/^ . 
We plot only points in the strong coupling region, b) Chiral equation 
of state on the 24^ lattice for the same Pc, S and Pmag-'m/ < ^''P >^ is 
plotted versus {Pc — P)/ < 4'^ ^'^/l^mag ^ y^jj ^]^g points are shown. 

7 Log-log plot of < ijjip > vs Mp. (Crosses, diamonds, circles, squares) are 

for p = (0.260, 0.255, 0.250, 0.245). 

8 < ip'^j ><^^-9l^ vsMp. pmag/v = (1,0.86/0.67,0.86/0.5) (dash, solid, dot 

line) 

9 Mp/{pc - PY vs. m/{pc - PY^ for v = 0.67, P = 0.84, 6 = 0.22. 

10 Mp{m = 0) vs. p. The sohd line is the best fit Mp = A{p - 0.260)°-^^^ 

11 logM^ vs. logm for p = (0.245,0.250,0.255,0.260,0.265) (top to bot- 

tom). The straight lines superimposed are fits to the relation M^ = 

12 /2 as a function of m ior p = (0.245, 0.250, 0.255, 0.260). 
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13 Scale invariant plot for {p/ fermionY vs (vr/p)^. (Squares, diamonds, 

crosses) are for /3 = (0.255,0.250,0.245). 

14 Scale invariant plot for [a / f ermionY vs (tt/p)^. (Squares, diamonds, 

crosses) are for p = (0.255,0.250,0.245). 

15 Scale invariant plot for (al/p)^ vs (tt/p)^. (Squares, diamonds, crosses) 

are for p = (0.255, 0.250, 0.245). 

16 Scale invariant plot for {cr/pY vs (vr/p)^. (Squares, diamonds, crosses) 

are for (3 = (0.255, 0.250, 0.245). 



26 



References 

[1] For recent reviews see for instance the contributions of J. B. Kogut and 
S. Love in Strong Coupling Gauge Theories and Beyond, ed. T. Muta 
and K. Yamawaki, World Scientific (1991) pp. 299 and 309 

[2] E. Dagotto, A. Kocic and J. B. Kogut, Phys. Rev. Lett. 60 (1988) 1250; 
61 (1988) 2416; 

[3] J. Bartholomew, J. B. Kogut, S. H. Shenker, J. Sloan, M. Stone, H. 
W. Wyld, J. Shigemitsu and D. K. Sinclair, Nucl. Phys. B230 (1984) 
222 

[4] V. A. Miransky, II Nuovo Cimento 90A, (1985) 189. P. Fomin, 
V.Gusynin, V.A. Miransky and Yu. Sitenko, Riv. Nuovo Cimento 6 
(1983) 1 

[5] C. N. Leung,S. T. Love and W. A. Bardeen, Nucl. Phys. B273 (1986) 
649; B323 (1989) 493. 

[6] W. A. Bardeen, S. Love and V. A. Miransky, Phys. Rev. D42 (1990) 
3514 

[7] A. Kocic , S. Hands,J. B. Kogut and E. Dagotto, Nucl. Phys. B347 
(1990) 217 

[8] B. Holdom, Phys.Lett B150 (1985) 301; R. Chivukula, A. Cohen and 
K. Lane, Nucl. Phys. B343 (1990) 554 and references therein. 

[9] V. Miransky, M. Tanabashi and K. Yamawaki, Phys. Lett. B221 (1989) 
177. 

[10] J .B .Kogut, E. Dagotto and A. Kocic, Nucl. Phys. B317 (1989) 253. 

[11] Y. Nambu and G. Jona-Lasinio, Phys. Rev 122 (1961) 345 

[12] A recent review on progresses in lattice QED is E. Dagotto in Dynamical 
Symmetry Breaking, ed K. Yamawaki, World Scientific (1992) p. 189. 

[13] E. Dagotto, J. B. Kogut and A. Kocic, Phys. Rev. D43 (1991) R1763 



27 



[14] A. Kocic,J. B. Kogut and M.-P.Lombardo, Chiral Symmetry, a — tt 
Splitting and PC AC, Illinois-CERN preprint, in preparation; A. Kocic, 
Phys. Lett. B281 (1992) 309. 

[15] D. J. Amit, Field theory, the renormahzation group and critical phe- 
nomena (Mc Graw-Hill, New York, 1978) 

[16] E. Dagotto, A. Kocic and J .B .Kogut , Nucl. Phys. B331 (1990) 500. 

[17] see for instance K.Bitar et. al, Nucl.Phys. B313 (1989) 348 

[18] K. C. Bowler et al, Nucl. Phys. B284 (1987) 299. 



28 



